'' ' ">~LL 



ACID PRECIPITATION IN ONTARIO STUDY 



LINEAR PROGRAMMING SCREENING MODEL 



FOR 



DEVELOPMENT AND EVALUATION OF 



ACID RAIN ABATEMENT STRATEGIES 



SEPTEMBER 1983 






TD 

195.54 

.06 

L56 

1983 




Ministry 
of the 
Environment 



! \J 



Copyright Provisions and Restrictions on Copying: 

This Ontario Ministry of the Environment work is protected by Crown 
copyright (unless otherwise indicated), which is held by the Queen's Printer 
for Ontario. It may be reproduced for non-commercial purposes if credit is 
given and Crown copyright is acknowledged. 

It may not be reproduced, in all or in part, part, for any commercial purpose 
except under a licence from the Queen's Printer for Ontario. 

For information on reproducing Government of Ontario works, please 
contact Service Ontario Publications at eopvright@oniario.ca 



ftftfrO 




•7?^/ 



LINEAR PROGRAMMING SCREENING MODEL 
FOR DEVELOPMENT AND EVALUATION OF 
ACID RAIN ABATEMENT STRATEGIES 



PREPARED FOR 



POLICY AND PLANNING BRANCH 
ONTARIO MINISTRY OF THE ENVIRONMENT 
135 ST. CLAIR AVENUE WEST 
TORONTO r ONTARIO r CANADA M4V 1P5 



BY 



Edward A. McBean and Associates Limited 

Waterloo, Ontario 

In Co-operation with 

Conestoga-Rovers & Associates Limited 

Waterloo, Ontario 



September 1983 
(Reprinted June, 1985) 



7929 



TABLE OF CONTENTS 



ENTERED JUL 1 5 2003 

MOE 

Standards Development Branch 
LIBRARY 



Page 

1.0 INTRODUCTION 1 

2.0 COST-EFFECTIVE APPROACH TO ABATEMENT OF 

ACIDIC DEPOSITION 3 

3.0 OVERVIEW OF MODEL STRUCTURE 8 

4.0 DATA INPUT REQUIREMENTS 11 

4.1 POLLUTANT SOURCES 11 

4.2 ABATEMENT TECHNOLOGIES 14 

4.3 RECEPTOR AREAS AND DEPOSITION TARGETS 20 

4.4 SOURCE/RECEPTOR POLLUTANT TRANSPORT 22 

5.0 MODEL OPERATION *ND OUTPUT 24 

5.1 INTRODUCTION 2 4 

5.2 LINEAR PROGRAMMING 24 

5.3 MODEL OUTPUT 2 7 

6.0 MODEL APPLICATION FOR STRATEGY ANALYSIS 29 

6.1 APPROACH TO THE ANALYSIS OF STRATEGIES 29 

6.2 EXAMPLE APPLICATIONS 31 

6.3 LIMITATIONS OF THE MODEL 3 4 

7.0 SUMMARY AND FURTHER APPLICATIONS 37 



APPENDIX I MATHEMATICAL MODEL DOCUMENTATION OF DATAGEN 
(under separate cover) 

APPENDIX II DEVELOPMENT OF SO2 EMISSION CONTROL 
COSTS 
(under separate cover) 

APPENDIX III CANADIAN SOURCE INVENTORY 
(under separate cover) 



1.0 INTRODUCTION 

The Ontario Ministry of the 
Environment has dedicated considerable time and resources to 
the study of the phenomenon commonly known as "Acid Rain". 
From this effort, impressive progress has been made in 
understanding the causes, effects and technologies available 
for reducing the emissions responsible for acid rain. 

This knowledge can provide the 
basis for determining acidic deposition objectives to 
protect the environment. It does not, however, solve the 
policy maker's problem of determining how much emission 
reduction, overall and by individual sources, must be 
undertaken to achieve specified deposition targets in 
receptor areas. Formulating and evaluating strategies for 
reducing acidic deposition is an arduous task, particularly 
because of the numerous approaches and criteria tnat may 
potentially be used to allocate emission reductions among 
SO2 emitters. For this reason, abatement strategy 
development is critically in need of a procedure to identify 
and "screen out" what are considered to be the most 
socially-desirable options. To this end, an analytical 
procedure was developed which reduces the number of possible 
policy options by identifying those strategies which are 
both efficient in terms of minimizing control costs and 
effective in terms of achieving the desired deposition 



objectives. Because the role of the procedure is to 
identify the most viable options, the model is referred to 
as the Screening Model . 

The purpose of this report is to 
provide a summary of the Screening Model. What follows is 
an overview; it is not intended that the overview be 
comprehensive in itself. Instead, a less formal discussion 
aimed at the general reader interested in the underlying 
logic of the approach is presented. For more detailed 
information describing the structure of the model, the 
reader is referred to the forthcoming appendices to this 
report (under separate cover) entitled, Appendix I "Mathe- 
matical Model Documentation of Datagen", Appendix II 
"Development of S0 2 Emission Control Costs", and Appendix III 
"Canadian Source Inventory". These reports will be available 
from: Policy and Planning Branch, 135 St. Clair Avenue West, 
Toronto. 



2.0 COST-EFFECTIVE APPROACH TO ABATEMENT OF ACIDIC DEPOSITION 

Traditional approaches to air 
pollution control were designed mainly to address local 
problems associated with high pollution concentrations near 
the pollution source. Achieving compliance witn ambient 
standards involved either tne source reducing its emissions 
or diluting the pollution concentration through 
construction of taller smokestacks. 

Unfortunately, existing air 
quality control programs are poorly suited to abatement of 
the acid rain problem. The focal point associated with acid 
rain damage is the total deposition loading over a region, 
as opposed to pollution concentrations. Moreover, acid 
pollution travels over great distances, often crossing 
jurisdictional boundaries, with the result that deposition 
at any location is the result of emission contributions from 
many sources. As a result of these distinctive 
characteristics, novel approaches for controlling acid 
pollution must be developed. Three general approacnes have 
been proposed: 

1) emissions-driven; 

2) deposition - targeted; and, 

3) cost-effective or least-cost targeted. 



An emissions-driven strategy refers to any program 
specifying emission reductions to be achieved by specific 
sources. Programs tnat would impose universal 50% emission 
reductions or alternatively reductions weighted by a 
sources' contribution to total provincial emissions, are 
examples of emission-driven abatement strategies. The major 
deficiency of such strategies is that source reductions are 
advocated independently of both their contribution to acidic 
deposition and their cost. 

The emissions-driven position, 
then, explicitly ignores fundamental aspects of acid rain, 
namely, that pollutants are transported over long distances 
and that acid deposition at any location is the result of 
emissions from many sources at different locations. 

Alternatively, a deposition- 
targeted approach considers the set of sources that must be 
controlled to achieve specified deposition levels at 
receptor locations. To accomplish this more informed 
approach to decision-making requires an understanding of the 
linkage or relationship, between source emissions and the 
resultant acidic deposition in defined receptor areas. Such 
linkage information is available through use of Long-Range 
Transport of Air Pollutant (LRTAP) Models that simulate tne 
complex physical and chemical processes which control 
pollution transport. Using LRTaP models, the contribution 



of any given emission source to the total acidic deposition 
received at any given receptor area can be estimated. 

From a policy context this 
capability is important on two accounts. First, it can be 
used to identify the acidic deposition consequences for 
receptor areas of implementing an emissions-driven strategy. 
Secondly, given a desired level of deposition for each 
receptor, the emission reductions necessary to achieve those 
objectives can be determined. 

Despite its improvements over the 
emissions-driven approach, the deposition-targeted method is 
still inadequate. Two problems are noteworthy. The set of 
emission source controls required to achieve specified 
deposition levels may not be unique. That is, a number of 
source-reduction combinations may yield the same deposition 
levels at receptor areas. None of these combinations may 
possess any inherently distinguishing feature from which a 
policy decision can be made. 

A further objection is that costs 
to control emissions are not considered. A regulatory 
agency which ignores control costs in tnis manner runs the 
risk of adopting policies which are not cost-effective. 



To redress these concerns, a 
methodology was developed which synthesizes cost of 
abatement and the source/receptor linkage information. This 
new orientation influences policy direction by providing 
answers to the questions: 

1 ) What effect does an abatement program have on deposition 
at receptor areas? 

2) What effect does the choice of strategy have on the 
costs of controlling emissions? 

Answers to these questions provide 
a basis for comparing the relative efficiency and 
effectiveness of abatement policies. As well, tne decision- 
maKers problem of having to select from a number of equally 
effective abatement options is eased by applying tne 
criterion of least-cost. The decision problem thus becomes 
one of searching for a set of emission reductions that will 
minimize abatement costs subject to the condition that the 
desired deposition levels are met. Under this least-cost- 
targeted approach, sources with the lowest marginal costs 
per unit reduction of deposition at critical receptor areas 
would undertake the greatest degree of control. 

The intent of this section was to 
briefly outline the rationale and advantages of adopting a 



least-cost targeted philosophy in establisning a policy for 
controlling acidic deposition. In the following section, 
the emission source inventory, source/receptor relationships 
and control costs which are input to the analytical model 
are discussed. 



3.0 OVERVIEW OF MODEL STRUCTURE 

Models are simplified 
characterizations of real systems tnat generally rely on 
mathematical expressions to describe and simulate processes 
and predict system behaviour. In the present model, three 
features help to capture the linkage between source 
emissions and receptor area deposition rates while enabling 
an analysis of cost-effective abatement strategies: 

i) Major sources and sensitive receptor areas (receptors) 
are individually identified and described; 

ii) Information from an atmospheric long-range transport 
model is used to link source emissions to receptor 
deposition levels; and, 

iii) Relying on the source-receptor linkage, a decision 

analysis technique called linear programming is used to 
identify levels of emission control at each source that 
will meet specified deposition targets at receptors 
while keeping emission control costs as low as 
possible . 

Within the modelled area, discrete 
point sources and diffuse sources generate SO2 emissions 
that, to a greater or lesser extent, contribute to the 



deposition of sulphates at each of the receptors. 
Long-range transport models* are used to simulate the 
physical and chemical processes that cause this 
contribution. Information from these long-range transport 
models, in the form of transfer coefficients**, is used in 
the Screening Model to link source emissions to receptor 
deposition levels. In effect, every source contributes to 
the deposition at each receptor, with the amount being 
explicitly accounted for in the Screening Model. This 
linkage is schematically depicted in Figure 1. 

The Screening Model is not a 
long-range transport model itself but, instead relies on 
information developed from such models. Rather, the 
Screening Model is designed expressly to analyze acid 
deposition management decisions. In this analysis, it 
relies on linear programming, a mathematical procedure that 
is used to "sort through" a large set of possible management 
options in order to identify one which is optimal from the 
point of view of specified criteria. These criteria are 
formulated either as constraints that must be met or as an 
objective to be achieved. For the acid deposition 



* e.g. Ontario Ministry of the Environment, "Documentation of 
the Ontario Ministry of the Environment Statistical LRT Model 
OME-LRT" , Work Group 2, Report 2-8, prepared for the 
Atmospheric Modelling Subgroup, March 31, 1981. 

** See Section 4.4 for a discussion of transfer coefficients. 
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FIGURE 1 - SOURCE/RECEPTOR 
LINKAGE IN THE SCREENING 
MODEL 



Screening Model, minimization of emission control costs can 
be formulated as the modelling objective. This objective is 
achieved subject to constraints requiring that receptor 
deposition levels be at, or below, maximum allowable 
limits . 

The emission control levels at 
each source are adjusted during model operation in order to 
satisfy both the objective and the contraints. They are the 
planning parameters about which decisions are being made 
within the model and thus are called decision variables. 

It is possible to include 
additional decision criteria in the acid deposition 
Screening Model such as might relate to policy or political 
considerations . 



10 



4.0 DATA INPUT REQUIREMENTS 



4.1 POLLUTANT SOURCES 



The data requirements for the 
model are substantial. The starting point for 
assembling the requisite data inputs is a network of 
sources and receptors. For the current model 
applications this network encompasses areas of Canada 
and the U.S. east of the Rocky Mountains. Emission 
sources within this area are introduced into the 
Screening Model either as large point sources or as 
aggregate "area" sources. 

A total of 235 large point sources 
have currently been selected on the basis that their 
1980 emissions exceeded 19,000 tonnes per year of 
SO2. Included in this group are Canadian and 
American power plants and non-power plants as follows: 

Power Plants Non-Power Plants Total 
Canadian 9 7 16 

American 178 41 219 



Total 187 48 235 

Estimated emissions for these 
large sources are: 
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Total 
Emissions 
for Sector 

(kilotonnes/yr) 

Large Power Plants 

in Canada 542.6 

Large Non-Power Plants 

in Canada 2169.0 



Percentage (%) of 
Country-wide 
Total* 



1 1 .8% 



47.3% 



Canadian Totals 4582.7 



59. 1% 



Large Power Plants 12721.7 

in U.S. 

Large Non-Power Plants 

in U.S. 1606. 1 



60.7% 



7.7% 



U.S. Totals 20954.8 



68.4% 



* Estimated for modelled areas in the respective 
country, i.e. for regions east of the Rocky 
Mountains . 



Overall, these large sources 
represent approximately 67% of total man-made emissions 
of SO2 that originated within the modelling area in 
1980. Detailed emission control data and associated 
cost information are currently given only for these 
large sources in the Screening Model. They are 
therefore the controllable sources. 

The remaining emission sources, 
consisting of a few large point sources, point sources 
emitting less than 19,000 tonnes per year, and diffuse 
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sources, have been aggregated by state and province 
into area sources. There are 92 such area sources, 
including : 

Power Plants Non-Power Plants Total Number 
Canadian 9 9 18 

American 32 42 74 

Total Number 41 51 92 

Area source emissions are 
currently fixed in that no control technologies are 
applied to them in the Screening Model. However, 
emission controls for area sources are being 
considered. They may be introduced to the model by 
assuming, for example, the desulf ur ization of neating 
oil at refineries and the resulting combustion of lower 
sulfur fuels at points throughout the area source. 

In addition to the emissions that 
originate from point and area sources, there are 
emissions both from natural sources and from man-made 
sources that lie outside of the boundaries of the study 
area. These sources are accounted for in the Screening 
Model by means of a "background deposition level" at 
each receptor. Currently, the background deposition 
levels have been set for all receptors at 2 kilograms 
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per hectare annually of wet-deposition sulfur, (6 kg 
S04/Ha-yr) but these values are individually assigned 
for the various receptors and therefore can be varied 
if necessary. 



4.2 ABATEMENT TECHNOLOGIES 

For each of the 235 large emission 
sources in the screening model, SO2 emission levels 
are specified along with estimated removal efficiencies 
and costs for new emission control techniques. 

Emission data are based on 1980 
estimates of emissions prepared by the Ontario Ministry 
of the Environment. These estimates reflect production 
activity levels for that year and the impact of control 
technologies that were in place at the time. 

Fuel sulfur content controls such 
as coal cleaning are handled at the plant level. For 
flue-gas and off-gas desulf ur ization techniques, 
removal efficiencies and control costs are determined 
for individual boiler and production lines that are 
found at each plant. Thus, a power plant with four 
coal fired boilers could nave limestone scrubbers 
assigned to each boiler separately as well as coal 
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cleaning for all boilers. The alternative emission 
control options at that plant would then consist of 
coal cleaning with or witnout the four scrubbers taken 
separately or in combination. Maximum reductions would 
occur with all four scrubbers installed plus coal 
cleaning . 

A certain amount of site-specific 
information was required to enable the analysis of 
SO2 emission controls. For power plants tnis 
includes : 

Item Units 

Annual Capacity Factor % 
Fuel - Type 

- Consumption Level weight 

- Sulfur Content % by weignt 

- Heat Content Btu/lb 
Boilers - Number 

- Capacities MW 

SO^ emission Tonnes/yr 

Existing SO2 Removal by Boiler % 

For each non-power plant emission 
source, information requirements were: 
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Item 

Operating Schedule 

Number of Emission Points 

Flue Gas Volume 

Flue Gas Temperature 

Flue Gas SO2 Content 

Production Process Type 

Existing S0 2 Removal by Point 



Units 
hours, days, weeks 

cubic feet/min 
°F 
% Volume 



Control technologies considered in 
the analysis of power plants include : 

i) physical coal cleaning (25% removal of coal 

sulfur ) ; 
ii) coal substitution down to 1.4% sulfur; and, 
iii) flue gas desulf urization using limestone scrubbing 

(up to 90% removal of SO2). 

For the industrial or non-power 
plant sources, a variety of technologies are 
considered, depending on the industry and production 
processes at each source: 



Process Name 



Limestone FGD 



Acid Plant 



Removal Application 

90% - weak gas streams from 
ooilers and various 
production processes 

95% - strong gas streams 
from non-ferrous 
smelters 
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wellman Lord FGU 90% - sulfuric acid plant 

off-gases 

Vacuum Carbonate 93% - cooking oven off-gas 

Process with Claus in the iron and steel 

Sulfur Recovery industry 

Petroleum Refinery Up to - petroleum refineries 

Controls 75% 



For each of the control techniques 
identified above, cost expressions were obtained that 
related variables such as fuel sulfur concentrations 
and flue gas flow rates to treatment cost. These 
expressions are based on published cost estimates for 
the various technologies. Separate expressions are 
used for capital costs and annual operating and 
maintenance costs. (See Appendix 2 for a detailed 
discussion of emission control costing.) 

The emission control cost 
expressions are in effect generic cost curves which 
fail to account for site specific conditions. They 
nave been applied to all sources for which site 
specific cost estimates were not available; namely all 
American sources and Canadian power plants outside of 
Ontario. For the remaining sources, more detailed cost 
estimates were obtained for S0 2 emission controls. 
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Total costs for a technology 
applied at an emission point in a plant were calculated 
as the sain of capital costs and the present value of 
annual operating and maintenance costs incurred over a 
20 year period. Capital costs represent the value at 
the end of the construction phase of expenditures 
incurred to put a control technology in place. 

The present value of future 
operating and maintenance expenditures is determined by 
adjusting future expenditures by an interest rate 
factor. This interest rate adjustment converts future 
expenditures into equivalent present day values thus 
accounting for the time value of money*. 

The present value of future 
expenditures measures the size of an annuity fund that 
would have to be created at the beginning of the 
planning period in order to just cover forecast future 
costs over the full planning horizon. The total cost 
therefore measures trie amount of money that would have 
to be provided in year one in order to cover 



* For an interest rate of 10%, the interest rate factor for a 
cost incurred in year "n" is: 
1/(1+. 10) n 

The present value of that cost is: P.V. = (cost in year n)/ 
(1 + . 10) n 



construction costs incurred up to that point in time 
and to cover all future operating and maintenance 
costs . 

Project costs in the Screening 
Model do not consider the possibility of production 
cost savings or by-product revenues. For the non-power 
plant sources these could at times be significant and 
an upward bias in certain costs could therefore occur. 
On the other hand, the generic cost estimates used here 
will be Low if unaccounted for factors act to increase 
costs . 

The primary analysis of removal 
efficiencies and costs associated with discrete 
emission points at a source provided a table of 
possible control options for that source as evidenced 
by the individual points plotted in Figure 2. From 
these control options taken individually and in 
combination, those which were most cost-effective, are 
identified as potential control measures to be input 
into the screening model. Berore being input to the 
screening model however, the plot of individual control 
options at a source was converted into a continuous or 
smooth removal-cost relationship. This is accomplished 
by joining key least-cost points with continuous linear 
segments as shown in Figure 2. 
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Legend: 
a, a -discretedata points 

▲ - least-cost points 

// - continuous curve 



50 
Percent Removal of SO2 



100 



FIGURE 2 - S0 2 REMOVAL/COST RELATIONSHIP 




FIGURE 3 - LOCATION 
OF SCREENING MODEL 
RECEPTORS 



The Screening Model uses the 
continuous removal-cost relationships in analyzing 
abatement strategies. In so doing, it is free to 
select any combination of cost and removal level found 
along the linear segments of the curve. Output from 
the Screeing Model is matched against the set of 
discrete removal-cost points to identify a specific 
control option that corresponds most closely to the 
point on the continuous relationship that was chosen by 
the model. (For further explanation, see Appendix I). 

Note that the concept of cost 
effectiveness is applied in deriving removal-cost 
relationships for each source as well as in the 
identification of an overall abatement strategy. At 
both levels of analysis, least-cost options are 
selected to control SO2 emissions. 



4.3 RECEPTOR AREAS AND DEPOSITION TARGETS 

Twenty locations are currently 
employed within the model to assess deposition rates in 
areas that are most susceptible to damage from acid 
deposition (see Figure 3). The model can be expanded 
to include additional receptors as needed. 
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The deposition criteria at these 
receptors are expressed in terms of maximum allowable 
receptor-deposition rates. These rates, which 
represent wet deposition only, are imposed as 
constraints upon management strategies generated by the 
Model. 

Currently, a single maximum 
allowable deposition rate is set for all of the twenty 
receptors, one particularly important level being 20 kg 
(S04)/ha-yr (or 6.7 kg S/ha-yr). This deposition is 
the sum of depositions attributed to: (1) each 
individual emission source, (2) each area source, and 
(3) the deposition caused by natural and unaccounted 
anthropogenic sources. Analyses to date indicate that 
a level of 20 kg (SC^/ha-yr) is only attainaole when 
emission controls are applied at a large number of both 
Canadian and American emission sources. Higher and 
lower target levels have been tested to check the 
response of emission control requirements to this 
parameter . 

One can also assign different 
target levels for deposition at each receptor if, for 
instance, assimilative capacities vary from region to 
region . 
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4. 4 SOURCE/RECEPTOR POLLUTANT TRANSPORT 

For all work to dare, the 
statistical model of the Ontario Ministry of the 
Environment has been used to produce transfer 
coefficients for the Screening Model. A transfer 
coefficient is needed for each source/receptor pair . 
As a result, 6,540 transfer coefficients are input into 
the model (i.e. 20 receptors x [235 point sources + 92 
area sources] ■ 6,540). 

Transfer coefficients represent 
unit response rates that measure the kilograms of 
sulfur deposited annually per hectare at the receptor 
that can be attributed to an annual emission of one 
tonne of SO2 at a source. 

The deposition at a receptor 
corresponding to tne base case emission from a 
particular source is estimated as: 

(base case emission in tonnes S02/year) x 
(transfer coefficient for that source-receptor 

pair) 

= wet deposition in kg S/ha-yr . 
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'Pile total wet deposition at a 

receptor is the sum of depositions from controllable 

sources plus deposition associated with area source 
emissions and background sources. 
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5 MODEL OPERATION AND OUTPUT 



5.1 INTRODUCTION 



The core of the Screening Model is 
the linear programming procedure. Computer programs 
were required to both prepare input data files to the 
linear programming procedure and to transform resulting 
output into a legible report. These two computer 
programming tasks entailed the development of a 
hierarchy of computer programs and data sets that 
together comprise the acid deposition Screening Model . 
This hierarchy is illustrated and briefly described in 
Figure 4. 

The primary data input files 
indicated in Figure 4 are the result of preparatory 
work in reviewing information sources and assembling or 
estimating data inputs. These files are processed by 
the input data generator programs to produce the input 
data sets for the linear programming procedure. 



5.2 LINEAR PROGRAMMING 

The linear programming decision 
variables are defined as the SO2 removal efficiencies 
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PRIMARY DATA INPUT FILES 








- 


background deposition levels 


by 


receptor 


— 


maximum allowable deposition 


by 


receptor 


- 


source/receptor transfer coe 


f f icients 


- 


point- and area- source 


emission 


data 


— 


removal- cost schedules 


by source 




INPUT DATA GENERATOR PROGRAMS 








- estimate uncontrollable deposition 


levels 


usi 


nq 


background depositions 


plus area source 






emissions 










- check feasibility of allowable depc 


sition 






rates 










- format data for linear 


programming 










LINEAR PROGRAMMING PROCEDURE 

- iterative search to identify the least-cost 
management strategy that satisfies allowable 
deposition constraints 




LINEAR PROGRAMMING SOLUTION 

- a minimum cost management strategy with 

emission control levels and associated costs 
given for each controllable source 




REPORT GENERATING PROGRAMS 

- estimate total deposition results by receptor 

- convert total costs into capital and annual 
operating/maintenance components 

- sum costs and removals by area 



FINAL SCREENING MODEL REPORT FOR ONE 
MANAGEMENT STRATEGY 

- deposition levels by receptor 

- removal requirements for SO2 emissions by 
state, province and country 

- removal costs by state, province and country 



FIGURE 4 - SCREENING MODEL STRUCTURE 
AND OPERATION 



at each controllable source (see Section 4.1 and 4.2). 
These removal efficiencies, expressed as a fraction 
between zero and one, enter into the equation 
describing the Model's least-cost objective where total 
removal costs are estimated. The coefficients in this 
equation measure removal costs. The removal efficiency 
variable multiplied times the cost coefficient will 
give the total emission control cost corresponding to a 
given level of S0 2 removal at a source. 

Total abatement costs are measured 
as the sum of removal costs incurred at every 
controllable source. This value is minimized using an 
iterative mathematical routine whereby all feasible 
combinations of control options selected aross the 
controllable sources are considered. There are a very 
large number of such combinations; however the 
iterative search using linear programming is relatively 
quick and efficient as a result of mathematical 
properties of the input data sets. 

In searching for the minimum value 
of the objective function, the linear programming 
routine progresses through two stages. First, costs 
are ignored and the Model merely adjusts the emission 
control removal levels until constraints such as those 
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describing allowable deposition rates are all 
satisfied. 

For example, the linear 
programming routine will increase the emission control 
removal efficiencies until total deposition rates at 
each receptor are no greater than the allowable rates. 
When this is achieved, the corresponding set of values 
for the emission control decision variables is called a 
"feasible solution". 

If the required reduction in 
deposition at any receptor exceeds the maximum 
attainable reduction that can be achieved then no 
feasible solution exists. In other words, the target 
deposition level cannot be achieved using the control 
technologies that underlie the removal-cost 
relationships in the Model. 

If a feasible solution can be 
found, the next stage in linear programming commences. 
This entails an iterative procedure to identify where, 
and by how much, emission control removal levels can be 
relaxed without violating deposition constraints. In 
this step, all of the many source/receptor interactions 
must be considered simultaneously. 



26 



Usually, there is only one 
least-cost feasible solution for any particular set of 
input data; this is called the "optimal solution". 



5.3 MODEL OUTPUT 

The optimal solution provides a 
data set prescribing emission control levels and 
associated costs for each of the 235 controllable 
sources. This data set is interpreted by report 
generating programs that produce final Model outputs. 

The primary Model outputs include 
the following: 

- removal efficiency (in percent) at individual 
controllable sources; 

- abatement costs for control at each source (capital, 
annual operating and maintenance, and present value 
of total costs in millions of U.S. 1981 dollars); 

- final emission levels at each source (tonnes SO2 
per year ) ; and, 

- deposition levels at each receptor arising after the 
individual source abatement efforts are in place (kg 
S per ha-yr) . Deposition levels will be at, or 
below, the target level. 
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For reporting purposes, several 
levels of output are available (see Figure 5). 
Detailed source data are aggregated by state, province 
and country to reduce the volume of output. Power 
plant and industrial data are reported separately. 
Deposition rate data for each receptor are also 
displayed. Examples of Model outputs are provided in 
Section 6.2. 
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FIGURE 5 -LEVELS OF AGGREGATION FOR SCREENING MODEL OUTPUTS 



6.0 MODEL APPLICATION FOR STRATEGY ANALYSIS 

6. 1 APPROACH TO THE ANALYSIS OF STRATEGIES 

The output from one application of 
the Screening Model provides a complete set of required 
abatement actions at each of the controllable emission 
sources, along with associated costs, and resulting 
deposition rates for a specified set of constraining 
conditions. This output comprises a data set for 
evaluating a single strategic option. 

Although there is only one 
"optimal solution" for any given set of input data, tne 
Screening Model is not meant to be used to identify the 
single most "optimal" management strategy. Even if 
such a strategy exists, it is unlikely to be discovered 
by any computer program in light of the data 
deficiencies and uncertainties that persist. Rather, 
the Screening Model functions to identify a small set 
of interesting or promising management strategies. 
This is achieved by successive applications of the 
Screening Model to input data sets that vary according 
to modifications in the underlying assumptions about 
the decision criteria or the behaviour of the acidic 
deposition system. 
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For example, deposition targets 
may be varied across receptors. This could be done to 
reflect varying assimilative capacities across 
receptors or to minimize the influence of certain 
receptors at which the deposition targets prove 
particularly difficult to meet. Alternative estimates 
of emission control costs, base case emissions or 
transfer coefficients may also be tested. 

Upper and lower bounds are 
specified for emission controls at each source. While 
these are typically used to reflect technological 
constraints on removal efficiency, they can be modified 
for instance to fix emissions for a source at existing 
levels (no additional control allowed) or to force a 
source into the management strategy (lower bound raised 
to prevent zero removal). In this way, complete sets 
of sources, designated by type or by political 
jurisdiction can be manipulated externally in order to 
explore different policy options. 

One can also use modifications in 
the input data set to look at the economic costs of 
imposing equity rules on the abatement efforts 
undertaken across political jurisdictions or across 
specified sectors of industry. One such rule would 
call for equal percentage reductions across nations . 
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By utilizing one or more of a wide 
array of possible constraints, the acid deposition 
Screening Model can be a very flexible tool in 
examining alternative management strategies. 



6.2 EXAMPLE APPLICATIONS 

Two model applications are 
described in the following pages. They are included 
only for purposes of explaining the type of results 
generated by the Model and are preliminary in nature 
since certain model inputs are still being updated. 
Both applications utilize the same "base case" 
structure or set of constraints; consisting only of the 
requirement that the maximum deposition rates are not 
exceeded at any of the twenty sensitive receptors. 
There are no "special" or additional constraints such 
as restrictions on removal levels within a political 
jurisdiction. The applications differ only in the 
maximum allowable deposition rate. Application number 
1 utilizes a target of 25 kg SC>4(wet) /ha-yr . The 
target for application number 2 is 20 kg 
SO4 (wet) /ha-yr . 

These two scenarios show how 
controllable source SO2 removal levels and 
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accompanying costs increase as the maximum allowable 
deposition target becomes more stringent (i.e. 
lower) . 

Before applying the Model, the 
maximum allowable deposition targets were modified to 
account for receptor-specific predicition errors 
inherent in the transfer coefficients obtained from the 
MOE Statistical Model. When compared to measured SO4 
deposition rates, the Statistical Model overpredicts 
deposition at some receptors while underpredicting it 
at others. These discrepencies can be expressed as 
"predicition factors" which are simply measured 
deposition rates divided by model-predicted deposition 
rates (deposition rates here being measured net of 
background deposition levels). The maximum allowable 
deposition rates used in the Model are the nominal 
values (25 and 20) divided by the receptor-specific 
prediction factors. 

Deposition results for the two 
applications are provided in Table 1. In accord with 
the fundamentals of linear programming, the search for 
minimum cost ensures that at least one receptor will 
have a deposition rate equal to the maximum. For a 
base target level of 25 kg SO4 (wet) /ha-yr, there 
are three such receptors: Muskoka, West Pennsylvannia 
and the Smokies. These receptors are termed "binding" 
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TABLE 1 

DEPOSITION RESULTS FROM TWO EXAMPLE APPLICATIONS 
OF THE SCREENING MODEL 



RECEPTOR 


NAME 


DEPOSITION VALUES 


(kg(S0 4 )/ha-yr) 


NUMBER * 


Target Level = 25 


Target Level = 20 


1 


Cree Lake, Sask. 


3.0 


2.9 


2 


Lynn Lake, Man. 


7.7 


7.0 


3 


Molson Lake, Man. 


8.7 


7.8 


4 


Goose Bay, Lab. 


10.4 


8.6 


5 


Red Indian Lake, Nfld. 


15.5 


13.7 


6 


Renous River , N . B . 


17.3 


14.3 


7 


Quebec North 


3.1 


3.0 


8 


Quebec Central 


8.2 


7.4 


9 


Quebec Southwest 


20.6 


16.3 


10 


Quebec South 


24.9 


19.9 


11 


Quebec Northeast 


17.6 


15.5 


12 


Boundary Waters 


6.9 


6.4 


13 


Algoma 


14.3 


11.0 


14 


Muskoka 


25.0 


18.3 


15 


Quebec City 


24.8 


20.0 


16 


Southern N.S. 


17.3 


13.8 


17 


New Hampshire 


24.2 


19.3 


18 


Adirondacks 


24.0 


17.7 


19 


West Penn. 


25.0 


16.2 


20 


Smokies 


25.0 


20.0 



* See Figure 3 for location of receptors. 



since the deposition targets here are the most 
difficult targets to satisfy. For this reason, the 
binding receptors largely govern the level of emission 
controls that are required. 

When the deposition target is 
reduced to 20 kg S04(wet) /ha-yr , the binding 
receptors become Quebec City and the Smokies. This 
change is evidence of the complexity of the 
source/receptor linkage; emission abatement at one 
source to attain maximum allowed deposition at one 
receptor also has an impact on deposition rates at 
other receptors. The ability to accommodate this 
complex interaction in an efficient manner is a 
principal value of the Screening Model. Analytic 
approaches that lack this capability must fall back on 
trial-and-error procedures that are inefficient and 
unreliable . 

The emission control levels and 
accompanying costs for the two applications are given 
in Table 2, for Canada and the United States. The 
aggregate removal levels and costs are higher for 
application 2 due to a more stringent deposition target 
in this application. 
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TABLE 2 

EMISSION CONTROL LEVELS AND COSTS FOR TWO EXAMPLE APPLICATIONS 

OF THE SCREENING MODEL 



1 



CANADA 



Power Non-Power 
Plants Plants 



UNITED STATES 

Power Non-Power 

Plants Plants Total 



Emissions (1Q3 Tonnes/yr) 
Existing (1980) 



542.6 



2169.0 



12721.1 1606.1 



17039.4 



After Controls 

a) Dep. target 

b) Dep. target 



25 kg S0 4 


419.1 


1384.1 


10451 .2 


1603.4 


13857.8 


20 kg SO4 


168.2 


644.2 


5407.5 


997.4 


7217.3 



Emission Reductions (%) 

a) Dep. target = 25 kg SO4 

b) Dep. target = 20 kg SO4 



22.7 
69.0 



36.2 
70.3 



17.8 
57.5 



0.2 
37.9 



19 
58 



Costs (1P6 U.S. 1981 Dollars) 

a) Dep. target - 25 kg S0 4 

Capital 

O&M 

Total 

b) Dep. target = 20 kg S0 4 

Capital 

O&M 

Total 



510.7 


581.5 


937.1 


0.4 


2029.7 


62.8 


50.6 


357.6 


0.1 


471.1 


1230.5 


1155.5 


4990.2 


1.6 


7377.8 






1844.3 


1025.6 


10773.2 


931 .8 


14574.9 


198.5 


132.4 


1698.7 


183.0 


2212.6 


4094.0 


2527.8 


30293.5 


3017.5 


39932.8 












The Screening Model output 
interpreter provides cost and removal results for three 
levels of aggregation; by country as shown in Table 2, 
by state and province, and individually by each 
controllable SO2 source. 



6.3 LIMITATIONS OF THE MODEL 

As with any modelling exercise, 
limitations of the Screening Model arise from sources 
of uncertainty and error in the input data and from 
constraints on the scope and complexity of the model 
structure. These limitations must be kept in mind when 
Model results are used to analyse management 
strategies . 

Invariably, data inputs are 
subject to a degree of uncertainty and error. For 
example, the base case emission levels from 1980 will 
clearly change as economic conditions shift and growth 
occurs. As scientific investigations proceed, data on 
transfer coefficients and deposition levels will 
improve. Such changes must have a bearing on our 
efforts to abate acid deposition. At this early stage 
in the strategy analysis, it is of paramount importance 
that sensitivity analysis be conducted in order to 
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assess the potential importance of possible future 
changes in the Model input data. 

All of the data inputs are subject 
to revision as better information becomes available. 
Indeed, the process of revision has already commenced, 
generic emission control costs for many Canadian 
sources have been replaced by more reliable 
site-specific emission control cost estimates. The 
degree of inference supported by Screening Model 
analysis must reflect the changing quality of the input 
data. While data inputs are in a very preliminary 
stage of refinement, Model outputs can be used to 
indicate relative rankings of management options and 
order of magnitudes of economic and environmental 
impacts. Absolute quantitative differences will take 
on more meaning only as the quality of the input data 
increases . 

Much the same can be said of Model 
scope and structure. It is easy to envisage an 
improved model encompassing additional pollutants 
(N0 X for example) , a more complete coverage of source 
and receptor areas, explicit measurements of receptor 
damages, and a more realistic scheduling of abatement 
activities. As it now stands however, many simplifying 
assumptions have been made. For instance, all 
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prescribed control options are assumed to come on line 
in year one. Consequently, no inferences about 
priorization of sources and scheduling of abatement 
activities can be made on the basis of Screening Model 
results alone. 

For the same reason, the financing 
and budgetory implications of abatement program costs 
cannot be assessed on the basis of Model results. At 
this stage, cost data in the model should only be used 
to evaluate the relative efficiency of alternative 
abatement strategies. 

Limitations cited above do not 
invalidate the Screening Model. The Model provides a 
powerful and practical means of synthesizing large 
amounts of complex and disparate information describing 
acid deposition. Existing limitations serve to define 
how Model outputs should be interpreted and they will 
help direct efforts for future Model development. 
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7.0 SUMMARY AND FURTHER APPLICATIONS 

The Screening Model described in 
this report is a computer model designed to assist in the 
discussion of strategies to abate acid deposition across 
eastern North America. In its operation it identifies a set 
of emission controls to be imposed at 235 major emission 
sources in Canada and the United States in order to achieve 
stipulated deposition targets at 20 receptors located in 
environmentally sensitive areas. The recommended controls 
represent the least cost manner of achieving the deposition 
targets. Only emissions of SO2 and wet depositions of 
sulfur are accounted for. A set of emission controls 
obtained from one application of the model provides the data 
with which to evaluate a single abatement program or 
management strategy. 

The principal merit of the model 
is that it provides, within a single model structure, the 
capability to simultaneously account for the complex 
interactions of multiple sources and receptors. Omission 
levels, control costs, pollutant transport and pollutant 
deposition at the receptors are all accounted for. The 
model therefore represents a significant advance in our 
efforts to develop an overview of the acid deposition 
problem. 
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Not only does the model foster an 
overview approach to the analysis of acid deposition, but it 
does so in a manner that is tailored to the needs of the 
decision maker. As a result of its reliance on linear 
programming/ a tool for decision analysis, it can address 
salient questions that bear on the management decisions; 
questions such as: 

- what can one political jurisdiction do on its own? 

- how far can abatement from power plants alone (or any 
other sector of the economy) take us? 

- how much more will it cost if everyone is forced to 
control emissions by say 50% rather than following the 
least cost option? 

Finally, the model is flexible. 
It does not have to be abandoned as new research findings 
and new data become available. Rather it can be readily 
modified and updated to incorporate new information. Indeed 
it is only as good as the data upon which it is based and 
must be revised if it is to remain a useful tool in the 
arsenal to combat acid deposition. 

This last point has implications 
for the future application of the Screening Model. Our 
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ability to comprehend the overall environmental system may 
diminish as the information available to us increases and 
becomes more complex. The Screening Model can help in this 
regard since the model structure can be adapted to 
accommodate new types of information. Two very pertinent 
examples of such information are receptor damages and 
transfer coefficient data for other acidic pollutants such 
as N0 X . Once estimates of receptor damages have been 
established for all of the significant receptor areas, then 
these can be entered directly into the economic calculations 
of the Screening Model. Qnission control costs could then 
be traded off directly against receptor damages in order to 
determine acceptable deposition targets. 

Apart from its role in the 
analysis of management strategies, the Screening Model may 
be able to provide some assistance in prioritizing acid 
deposition research efforts. Much of the information on 
acid deposition that is currently available is subject to 
varying degrees of uncertainty. Using sensitive analysis, 
the significance of these various sources of uncertainty to 
decisions regarding abatement can be assessed with the 
Screening Model. Research efforts can be directed to those 
areas where uncertainty appears to be most critical. 

Finally, it is conceivable that in 
a revised form, with significantly improved input data, the 
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Screening Model could be used daring the implementation 
phase of an acid deposition abatement program. It could 
provide one means of assessing, for example, the economic 
implications of various degrees of control at important new 
sources or of evaluating economic penalties for 
non-compliance on the part of polluters. 
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